RELAXED STATES OF 
A MAGNETIZED PLASMA 
WITH MINIMUM DISSIPATION 



B. Dasgupta*, P. Dasgupta**, 
M. S. Janaki*, T. Watanabe*** and T. Sato*** 
*Saha Institute of Nuclear Physics 
I/AF, Bidhannagar, Calcutta 700 064, India 
**Department of Physics, University of Kalyani, Kalyani 741235, India 
* * * Theory and Computer Simulation Center 
National Institute for Fusion Science, Toki, Gifu 509-52, Japan 



Abstract. 



Relaxed state of a slightly resistive and turbulent magnetized plasma is 
obtained by invoking the Principle of Minimum Dissipation which leads to 

VxVxVxB = AB. 

A solution of the above equation is accomplished using the analytic continu- 
ation of the Chandrasekhar-Kendall eigenfunctions in the complex domain. 
The new features of this theory is to show (i) a single fluid can relax to an 
MHD equilibrium which can support pressure gradient even without a long- 
term coupling between mechanical flow and magnetic field (ii) field reversal 
(RFP) in states that are not force-free. 
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In the well-known theory of relaxation of magnetoplasma, Taylor |1| pro- 
posed that the process of relaxation is governed by the principle of mini- 
mum total magnetic energy and invariance of total (global) magnetic helic- 
ity K — J v A ■ BdV where the integration is over the entire volume, the 
latter being the most significant invariant in the theory of relaxation. Ac- 
cordingly the relaxed state of a magnetoplasma satisfies the corresponding 
Euler-Lagrange equation 



V x B = AB (1) 

with constant A, and, consequently, is a force-free state. Taylor's theory is 
quite successful in explaining a number of experimental results, including 
those of RFP. However, relaxed states as envisaged by Taylor, have only zero 
pressure gradient. 

Extensive numerical works by Sato and his collaborators have estab- 
lished @,@ the existence of self-organized states with finite pressure, i.e. 
these states are governed by the magnetohydrodynamic force balance rela- 
tion, namely, j x B = Vp, rather than j x B = 0. Recently, it has been 
demonstrated both by numerical simulation J|] and by experiments [f| that 
counter-helicity merging of two spheromaks can produce a Field- Reversed 
Configuration (FRC). The FRC has zero toroidal magnetic field and the 
plasma is confined entirely by poloidal magnetic field. It has a finite pressure 
with a relatively high value of j3. It may be concluded that FRC, with its 
non zero perpendicular component of current, is a relaxed state and it is a 
distinctly non-force free state. From the point of view of plasma relaxation, 
the formation of FRC through the counter-helicity merging of two sphero- 
maks is a unique process where a non-force free state emerges from the fusion 
of two Taylor states. The conclusion is that there exists a general class of 
relaxed states which are not always force-free, and Taylor's force-free states 
constitute only a subclass of this wider class. While Taylor states do not 
support any pressure gradient, equilibrium obtained from the principle of 
minimum energy accommodates pressure gradients only in presence of flow. 
Several attempts||,||] have been made in the past to obtain relaxed states 
which could support finite pressure gradient, a large number of them making 
use of the coupling of the flow with magnetic field|5|]-JT2]]. 

The principle of " minimum rate of entropy production" , formulated by 
Prigogine |T3[ and others, is believed to play a major role in many problems of 



irreversible thermodynamics. Dissipation, along with nonlinearity, is ubiqui- 
tous in systems which evolve towards self-organized states. Another closely 
related concept, the principle of minimum dissipation rate was used for the 
first time by Montgomery and Phillips |14| in an MHD problem to understand 



the steady state profiles of RFP configuration under the constraint of con- 
stant rate of supply and dissipation of helicity the usual physical boundary 
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conditions for a conducting wall. It may be pointed out that the principle of 
minimum dissipation was also discussed by Chandrasekhar and Woltzer [|16| 



in a sequel to the complete general solution of the force-free equation by 



Chandrasekhar and Kendall [15]. The minimum dissipation rate hypothesis 



was later used by a number of authors []17l| , fl8fl to predict the current and 
magnetic field profiles of driven dissipative systems. 

This paper deals with the question of determining the field configurations 
assumed by a magnetofluid in a relaxed state in absence of any external 
fields, while maintaining that the relaxation is governed by the hypothesis of 
minimum rate of energy dissipation. It is our conjecture that relaxed states 
could be characterized as the states of minimum dissipation rather than 
states of minimum energy. The novel feature of our work is to show that it 
is possible for a single fluid to relax to an MHD equilibrium with a magnetic 
field configuration which can support pressure gradient, even without a long- 
term coupling between the flow and the magnetic field. In a recent work, 
Steinhauer||12|| has claimed that single fluid MHD theory can admit only a 
force free state and one need to take recourse to a two fluid theory so that 
electromechanical coupling produces pressure gradient and a non force free 
state. Our work establishes that none of these requirements need be satisfied 
to obtain a relaxed state of the desired kind. 

In what follows we derive the Euler-Lagrange equation from a varia- 
tional principle with minimum energy dissipation and conservation of total 
magnetic helicity, solve the equation in terms of the analytically continued 
Chandrasekhar- Kendall eigenfunctions, discuss the important role played by 
the boundary conditions, and present our results for the flux, field reversal 
parameter and pinch parameter. We also compute the helicity integral, and 
show the plots of magnetic field, current, and pressure profiles. The field 
reversal parameter from our theory is definitely in better agreement with the 
experimental value than what is obtained from Taylor's theory. 

We consider a closed system of an incompressible, resistive magnetofluid, 
without any mean flow velocity, described by the standard MHD equations 
in presence of a small but finite resistivity n. In the absence of any externally 
imposed electric fields, the ohmic dissipation rate R is itself a time varying 
quantity. However, it is possible to find constraints that are better preserved 
than the rate of energy dissipation, so that the system self-organizes to certain 
relaxed states, which remain stable over time scales short compared to ohmic 
dissipation time. In this case, helicity still serves to hold as a good constraint 
as it decays at a time scale much slower in comparison to the decay time scale 
of the rate of energy dissipation as is evident from the simulation works of 
Zhu et. al. Q. In the following, we compare the decay rates of the energy 
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dissipation rate and helicity. The former is obtained as 



— - 9?f- V k A h 4 

dt s 2 v 

while the latter turns out to be 

dt S k 

From the above two equations, we see that the decay rate of energy dissi- 
pation is once again 0(1) at scale lengths for which k ~ S%. But at these 
scale lengths, helicity dissipation is only 0(S' _1 / 2 ) << 1. Thus, we may ex- 
pect that in presence of small scale turbulence, the rate of energy dissipation 
decays at a faster rate than helicity. 

We therefore minimize the ohmic dissipation R = J r]j 2 dV subject to the 
constraints of helicity / A • ~QdV . The variational equation is given by 

r]f + AA • B) dV = (2) 

where A is Lagrange's undetermined multiplier. The variation can be shown 
to lead to the Euler-Lagrange equation 

V x V x V x B = AB (3) 

where, A = A/?7 is a constant. The surface terms in the equation vanish if 
we consider the boundary condition 5 A x n = as well as j x n = 0, which is 
the physical boundary condition we will impose in the problem. 

We like to emphasize that eq. (3) is a general equation which embraces the 
Woltzer- Taylor equation ( i.e. eq. (1) ) as a special case. Now we proceed to 
construct a solution of eq.(3) and show that the general solutions can have 
j x B / 0. In other words, eq.(3) may lead to a non force-free state. 

The solution of eq.(3) can be constructed using the Chandrasekhar- Kendall 



(CK) eigenfunctions. Chandrasekhar and Kendall's solution ]n| of the equa- 
tion V x B = A B can be written (with three parameters //, m, k in cylindrical 
coordinates), as, 



B(/x, m, k) = AV$ xVz + Vx (V$ x Vz) (4) 

where, $ = J m (/ir) exp[i(m9 — kz)) with A 2 = /i 2 + k 2 . Here, J m is a Bessel 
Function of order m and the value of /i in the argument is determined from 
the boundary condition at r = a, which is given as (n • B) r=a = 

Analytic continuation of the above solution for complex values of \x (or k) 
is straightforward. For real values of A (and hence of \x and k) the operator 
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(Vx) has been proved to be self-adjoint, but not so in the larger space 
spanned by the analytically continued CK solutions. 
We introduce the complex parameters 

r 1 1/2 

fx n = (y + k 2 )exp(An-K i/3) - k 2 , n = 1, 2 (5) 
so that fi 2 +k 2 = \ 2 u 2n , uj = exp(2ni/3), and define 

Bi = B(/i, m, k) = AV$ xVz + Vx (V$ x Vz) 

B 2 = B(//i,m, fc) = AcjV$i x Vz + V x (V$i x V*) 

B 3 = B(/i 2 , m, fc) = Aw 2 V$ 2 x V2 + V x (V$ 2 x Vz) (6) 

In the last two expressions above, $i and $ 2 are obtained from $ by replacing 
fi by fii and //2 respectively. 

A solution of eq.(3) can now be obtained as a linear combination of 
Bi, B 2 , B 3 : 

B = «iBi + a 2 B 2 + a 3 B 3 (7) 

where a, are constants, with at least two of them non-zero. It can be easily 
demonstrated that the expression for B given in (7) is a solution of eq.(3) 
with A = A 3 . 

A reasonable boundary condition is to assume a perfectly conducting wall, 
so that 

B • n = 0, j x n = at r = a (8) 

The boundary conditions given by eq. (8) suffice to fix the arbitrary 
constants 

(9) 
(10) 

The magnetic fields at the boundary r = a have to obey the following 
relation for non-trivial values of the constants a, 

2B lr Im(B 2e B* 2z ) - 2B ie Im(iu 2 B 2r B; z ) + 2B lz Im(u 2 B 2r B* 29 ) = (11) 

From eq. (6), it is evident that B 2 and B 3 are complex conjugate of one 
another. This, together with the relations obtained in eq. (10), shows that 
the magnetic field given by eq. (7) is a real field. We also list the follow- 
ing expressions for the m = 0, k = state (cylindrically symmetric state) 



a 2 


{jj 2 {B xe B* 2z 


B 2 qB\ z ) \ r=a 


(X\ 


{B 2e B^ z - 


B 2e B 2z ) \ r=a 


a 3 


„ * 

= OL 2 
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obtained from eqs. (4)-(7) : 



B r 

Br 





A 2 «i 



Ji(Ar) + 2Re (^V Ji(Xojr)^ 
J (Ar) + 2Re (^uj 2 J {Xujr^j 



For a given value of m and fca, the value of Xa can be obtained from 
the boundary condition given by eq. (11). It is to be noted that for the 
cylindrically symmetric state the boundary condition is trivially satisfied 
and hence does not determine Xa. It can be easily proved that the state of 
minimum dissipation is equivalent to the state of minimum value of A. To 
get the numerical value of A for m ^ 0, we solve numerically eq.(ll) and 
obtain Xa = 3.11 and ka — 1.23 as the minimum values of Xa and ka for the 
m — 1 state. 

The only undetermined constant in eq. (7) is the value of a.\ (the value of 
the field amplitude) which can be determined by specifying the toroidal flux 
$ z . The m = k = state is responsible for non-zero values of toroidal flux 
which is obtained as 



$ 2 = 2naiXa 



Oil 

JAXa) + 2Re\—ujJ 1 {Xuja)} 



(12) 



A couple of dimensionless quantities that have proved useful in describing 
laboratory experiments are the field reversal parameter F = B z (a)/ < B z > 
and the pinch parameter = Bg(a)/< B z >, where < .. > represents a 
volume average. After substituting the expressions for B z (a) etc, we get 



F 




Xa J (Xa) + 2Re [(a 2 / ai)uo 2 J (Xuja)] 
2 Ji(Aa) + 2Re [(a 2 /ai Wi(Aa>a)] 

Xa Ji(Xa) + 2Re [(a 2 / «iV 2 Ji(Xua)] 
2 Jx(Aa) + 2Re [(a 2 /a 1 )ujJi(Xuja)] 



(13) 
(14) 



The pinch ratio is related to the ratio of the current and flux and is a 
physically controllable quantity. For the Taylor state = Aa/2. 

The details of any relaxed state are determined by two physically mean- 
ingful parameters, the toroidal flux and the volts-seconds of the discharge. 
The toroidal flux as defined earlier serves to determine the field amplitude 
and the volts-seconds describes the helicity of the relaxed state through the 
relation: volts — sec = helicity /flux 2 . We therefore calculate the helicity in- 
tegral (global helicity) from our solution for the m — 0, k = state and 
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get 



K = 



An 2 al^{\af J 2 {Xa) + J 2 (Xa) 

2 

2Re ^| Jl{\uoa) + J^(Xujo) - 



2 

Xa 



Jo(Aa) Ji(Aa) 



2 



J (Xua) Ji(Xuja) 



+ 



Xuoa 



+ 



+ 




ujJ (Xtu 2 a)J 1 (Xuja) — J 1 (XLU 2 a)J (Xuja) 



[Jo(Xa) Ji(Xcua) — uJ\(Xa) Jo(Xua)] 



(15) 



We then calculate the volts-seconds for the cylindrically symmetric state 
using Xa = 3.11. For the minimum value of Xa = 3.11, the critical value of 
volts-seconds = 12.8 R/a. For values of volts-seconds less than this critical 
value, a lower value of Xa is obtained from solving the equation for K/<& 2 
so that the cylindically symmetric state is the relaxed state for minimum 
energy dissipation. For values of volts-seconds greater than the critical value 
the system relaxes to the helically distorted state with Xa = 3.11 which is 
obtained as a mixture of the m = 0, k = and the m — 1, k ^ states as in 
the case of Taylor's theory. 

The profiles for the current and magnetic field are shown in Figs. 1 and 
2 for the m = 0, k = state with the value of Xa < 3.11 (i.e., volts-sec 
< 12.8i2/a). At a value of Xa greater than 2.95, the magnetic field profile B z 
vs r shows a reversal near the edge ( Fig. 2 ). Also, j z and jg go to zero at 
the wall because of the boundary conditions we have chosen. The values of 
both F and © at the boundary r = a are evaluated and F is plotted against 
pinch ratio ( Fig.3 ). It is observed that F reverses at a value of = 2.4, 
(Xa = 2.95) whereas for the Taylor state the reversal is achieved at = 1.2. 
However, this field reversed state supports pressure gradient in constrast to 
the Taylor state. The q-profile, where q = rB z /RB e is shown in Fig. 4 for 



The pressure profile can be obtained from the relation j x B = Vp. For 
the m = 0, k = state, the only nonvanishing component of the pressure 
gradient exists in the radial direction. The pressure profile is shown in Fig. 
5 for the m = k = state with Xa = 3.0 which is the minimum energy 
dissipation, field reversed state. 

To conclude, the principle of minimum dissipation is utilized together 
with the constraints of constant magnetic helicity to determine the relaxed 
states of a magnetoplasma not driven externally. The variational principle 
leads to a remarkable Euler-Lagrange equation, and it is shown that this 
equation involving higher order curl operator can be solved in terms of of an 
analytical continuation of Chandrasekhar-Kendall functions in the complex 



Xa = 3.0. 
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domain with appropriate boundary conditions. This relaxed state obtained 
from single fluid MHD supports pressure gradient. A coupling between mag- 
netic field and flow is not an essential criterion for having a non-zero pressure 
gradient. Further, it is shown that a non force-free state with field reversal 
properties can exist. 
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Figure Captions 



Fig.l jzij') vs r for the axisymmetric state m — 0, k — 0, and Aa = 3.0. 
The current vanishes at the edge because of the boundary condition j x n = 0. 

Fig. 2 Magnetic field profile for the axisymmetric state with Aa = 3.0, 
showing field reversal near the edge. 

Fig. 3 The field reversal parameter F against the pinch parameter G, 
the field reversal occuring at = 2.4. The dotted curve represents the plot 
for the minimum energy state of Taylor. 

Fig. 4 The q-profile for the axisymmetric state. 

Fig. 5 The pressure profile p vs r for Aa = 3.0. 
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